Geometric local invariants and pure three-qubit states 
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We explore a geometric approach to generating local 5(7(2) and SL{2, C) invariants for a collection 
of qubits inspired by lattice gauge theory. Each local invariant or 'gauge' invariant is associated to 
a distinct closed path (or plaquette) joining some or all of the qubits. In lattice gauge theory, the 
lattice points are the discrete space-time points, the transformations between the points of the lattice 
are defined by parallel transporters and the gauge invariant observable associated to a particular 
closed path is given by the Wilson loop. In our approach the points of the lattice are qubits, the 
link-transformations between the qubits are defined by the correlations between them and the gauge 
invariant observable, the local invariants associated to a particular closed path are also given by 
a Wilson loop-like construction. The link transformations share many of the properties of parallel 
transporters although they are not undone when one retraces one's steps through the lattice. This 
feature is used to generate many of the invariants. We consider a pure three qubit state as a test 
case and find we can generate a complete set of algebraically independent local invariants in this 
way, however the framework given here is applicable to generating local unitary invariants for mixed 
states composed of any number of d level quantum systems. We give an operational interpretation 
of these invariants in terms of observables. 

PACS numbers: 03.67.Mn 



I. INTRODUCTION 

One approach to the study of entanglement is the iden- 
tification of local invariants of a collection of quantum 
objects. With this approach we imagine the distant labs 
scenario in which N spatially separated parties each hold 
one of the subsystems of an N particle entangled state 
in their laboratory and they are free to make arbitrary 
transformations on their subsystem. One then looks for 
properties of the state that remain unchanged under such 
local transformations since, under the conditions that 
the transformations are unitary, entanglement is defined 
to be invariant. If the transformations belong to the 
group SL{2, C) it turns out that entanglement, given by 
the well known measure concurrence, is also invariant. 
Rephrasing, we can write this scenario as a non-abelian 
lattice gauge theory; the arbitrary transformations are 
non-abelian local gauge transformations made on N sub- 
systems, the N points of the lattice. Entanglement is a 
gauge invariant observable of the theory. It is this simi- 
larity that inspires our work. 

Quite a lot is known about the local unitary invari- 
ants of simple entangled states. For example, for a pure 
state of a pair of qubits, there is essentially only one local 
invariant (not counting the normalization); it character- 
izes the amount of entanglement between the two qubits. 
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For a pure state of three qubits, one can identify five 
independent local invariants, four of them expressing a 
different aspect of the state's entanglement [ll-H. A fifth, 
the Kenipe invariant, is not well understood [4']. There 
exist well known algebraic methods for generating invari- 
ants but as the number of subsystems increases, 
the problem of identifying and interpreting independent 
invariants rapidly becomes very complicated. 

Here we explore a different approach inspired by lat- 
tice gauge theory For a collection of N qubits, we 
consider any closed path connecting some of the qubits, 
and we associate an invariant quantity with each such 
path. The invariant is formed by taking the trace of a 
transformation associated with the closed path, which in 
turn is the product of transformations associated with 
the individual two qubit links. Each of these 'link- 
transformations' is determined by the density matrix 
of the two qubits connected by the link. Because this 
density matrix will typically change if one performs a 
local operation on either of the two qubits, each link- 
transformation will also typically change under such local 
operations. The overall transformation around a closed 
loop can similarly change as one performs a local opera- 
tion on the qubit that defines the loop's starting point. 
However, the trace and the eigenvalues of the overall 
transformation do not change under any single-qubit op- 
erations. The trace is our invariant. In fact we will gen- 
erate a few distinct invariants associated with the same 
path, by using different, but closely related, ways of mak- 
ing the correspondence between a two-qubit density ma- 
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trix and a link-transformation i.e. one has the choice 
whether to apply a spin flip to each qubit in a given 
loop. 

Other authors have explored relations between entan- 
glement and gauge transformations, in the context of an 
analysis of the geometry of the set of states [l^ - [l5j | . Our 
approach is different in that the paths we consider are not 
paths in the set of states but discrete paths connecting 
the qubits themselves. 

Thus our invariants are determined once we specify the 
correspondence between a two qubit density matrix and 
a link-transformation. The first rule we consider, and the 
one from which the other cases are derived, is this one: 



Mb = tra [{Ma <E) Ib)pab] 



(1) 



Here pab is the two qubit density matrix in question. Ma 
is a 2 X 2 Hermitian matrix, and Mb is its image under 
our transformation. (In section we interpret this rule 
in terms of local observables.) 

We hope that this geometric approach will ultimately 
prove useful in generating and classifying invariants of 
systems with many parts. In this paper we try out our 
ideas by applying them to a simple system of three qubits 
in a pure state. For that indicated above, a 

natural set of SU{2) invariants is already known. We ask 
whether this set, or an equivalent set, can be generated 
via our construction. We also ask whether the path-based 
approach sheds any light on the physical meaning of these 
invariants. 

In the following section we introduce our basic path- 
based method of generating invariants. Section IIVI ap- 
plies this idea to pure states of three qubits and makes 
the connection between the invariants generated by this 
method and the standard three qubit SU{2) invariants 
that have been identified previously. In section |V] we 
show how to generate SL{2, C) invariants by simply spin 
flipping each qubit in a loop. In section I VII we give an 
operational interpretation of these invariants in terms of 
observables. Finally, we draw conclusions in section IVIII 
and outline how one would extend this approach to mixed 
states comprised of any number of qudits. 



II. PATH-BASED INVARIANTS 

Our basic method of associating a transformation with 
each two qubit link is motivated by a thought experiment. 
Imagine many N qubit systems, each having distinguish- 
able qubits labeled a, 6, c, . . ., and each system being in 
the same quantum state p. We use p to define a trans- 
formation from qubit a to qubit h as follows. On several 
copies of the state p, perform a general quantum measure- 
ment on qubit a, one of whose outcomes is represented 
by the operator Ma- (This operator is arbitrary except 
that it must be positive semi-definite and less than the 
identity so that it can be part of a legitimate measure- 
ment.) Now consider only those instances of qubit a for 
which this particular outcome is actually achieved. In 



those cases, the state of qubit h has been 'collapsed' into 
some state, typically mixed, even though qubit h has not 
interacted with the measuring device. The final state of 
qubit b is in fact proportional to 



Mb^ira[{Ma®lb)Pab] 



(2) 



where pab is the original reduced density matrix of qubits 
a and h when the whole system is in state p. The normal- 
ization of Mb, that is, trM^, is equal to the •probability of 
getting the outcome represented by Ma- In this way the 
density matrix pab defines a linear transformation from 
operators on qubit a to operators on qubit 6, namely, the 
transformation that takes Ma to Mb- It is convenient to 
represent this linear transformation as a matrix by writ- 
ing Ma and Mb in terms of Pauli spin matrices. Let the 

0, . . . , 3, be defined by 



four real numbers m", i 



Ma = mgcTo -I- mlai + m2(J2 + ra^a^i, 



(3) 



where ctq is the 2x2 identity matrix and the other UiS 
are the usual Pauli matrices, and let the components of 
Mb be defined similarly. Then we can express our trans- 
formation as the 4x4 matrix S{b, a) such that 



m'' = S{b,a)m'', 



(4) 



where m° and m'' are column 4- vectors with components 
to" and m^. Writing out eq. Q explicitly in this operator 
basis, we have 



^ TO^CT^ = tra 



Pab 



(5) 



Multiplying both sides by tj^ and tracing over qubit 6, 
we get an explicit expression for the components of the 
matrix S{b, a)ji 



S{b,a) 



1 

= -tr I 
2 



cr'^jPab) 



(6) 



So in this representation, the matrix representing our 
transformation from qubit a to qubit b is proportional to 
the spin correlation matrix. Our link-transformations are 
specified by the correlations between the qubits joined by 
the link. 

We can now imagine repeating this process at qubit b. 
That is, starting with several pristine, unmeasured copies 
of the state p, we imagine performing a measurement on 
qubit b, one of whose outcomes is represented by the same 
operator Mb that was the result of the first measurement. 
When this outcome is achieved, qubit c will be collapsed 
into some mixed state proportional to Mc defined as in 
eq. @. 

Continuing in this way around a closed loop, we finally 
collapse qubit a into some state proportional to 



= tr, [{M, (g) Ia)Pza] 



(7) 



where qubit z is the one that precedes qubit a at the end 
of the loop. In this way we have mapped, via the whole 



3 



S(b,a) 



-> — • — 



FIG. 1: A schematic of our method of generating local invari- 
ants or 'gauge invariants' for a collection of quantum objects. 
The idea is inspired by lattice gauge theory with the lattice 
points representing qubits and the transformations S{h, a) be- 
tween lattice points a and b given by the correlations between 
qubits a and h. Each distinct closed path or loop generates a 
local invariant. 



loop C, an operator Ma on the state space of qubit a into 
another operator M'^ on the same space. The matrix 
representing this transformation is 



S{a,a]C) — S{a, z) ■ ■ ■ S{c, b)S{b, a). 



(8) 



In other words the overall transformation taking our ini- 
tial measurement four vector m" around the closed loop 
back to our new four vector m° is 



m° = S'(a,a;C)m^ 



(9) 



As we show in the following section, the trace of this 
matrix is invariant under all single qubit unitary trans- 
formations. This trace is the invariant we associate with 
the given closed path. We present a graphical illustration 
of the idea in figure [TJ 

Our basic transformation, eq. ([2]), is reminiscent of the 
transformation that would be associated with the two 
qubit state pab by the Jamiolkowski isomorphism [l6l |. 
which provides a general correspondence between bipar- 
tite states and trace-preserving operations. The trans- 
formation defined by that isomorphism would be 



2tr a[{M^ (g) lb) Pab] 



(10) 



That is, it would be normalized differently and it would 
require taking the transpose of the initial operator. The 
transpose is included in order to make the transformation 
a legitimate quantum operation — specifically, in order to 
make it completely positive. In contrast, the transfor- 
mation defined by our eq. ([2]) need not be completely 
positive. We have chosen the form of eq. ([2|) as we have 
because our aim is not to define a quantum operation 
but rather to generate invariants. If we had included 
the transpose, the resulting trS'(C) would not have been 
invariant. Moreover, even though our transformation is 
not a quantum operation, it does have a physical inter- 
pretation in terms of measurement. 

Though introducing the transpose would spoil the in- 
variance, there is a closely related operation that does not 
have this effect, namely, the spin flip. At any point along 
a closed path, we have the option of inserting a spin flip 



without destroying the invariance. In our thought ex- 
periment, this would mean that, having collapsed, say, 
qubit b into the (subnormalized) state Mb, in our next 
step we would perform a measurement with an outcome 
represented by Mb, where the tilde represents the spin 
flip. That is. 



M = a2M^ (72- 



(11) 



The effect of a spin flip on the vector m representing M 
in the Pauli basis is simply to multiply mi, m2, and 7713 
by —1 and to leave toq unchanged. ((T2 multiplies mi and 
TT^3 by —1 while the transposition, equivalent to complex 
conjugation, multiplies m2 by —1. The spin flip is anti- 
unitary and not a physical operation). That is, in this 
representation a spin flip is represented by the matrix rj, 
the Minkowski metric; 



(12) 



We will label each of our invariants by the closed path 
that defines it, indicating with a tilde any site at which 
we have added a spin flip. Thus, for example, I{abc) is 
the invariant deflned by 
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I{abc) = tr {S{a, c)S{c, b)r]S{b, a)} . 



(13) 



It is not hard to see that a spin flip indeed preserves the 
SU{2) invariance. In the following section, we will show 
local transformations U become elements of S0{3) acting 
only on the spatial dimensions of S, those with index 
values 1,2,3 and not on the dimension associated with 
the identity. That is, they are block-diagonal matrices 
with a 1 X 1 block and a 3 x 3 block. Thus they commute 
with 77 and therefore still cancel each other. 

In fact we find that the inclusion of a spin flip on every 
qubit in a path results in not only an SU{2) invariant but 
also a SL{2, C) invariant, a group which contains SU{2). 
The group SL{2, C) represents the most general, local 
operations, such as Kraus operations, that one may per- 
form on a qubit up to a positive constant less than unity. 
This stronger invariance is interesting as the well known 
entanglement measures concurrence and three-tangle ex- 
hibit this higher invariance [l7j . We demonstrate this in 
the following section. 



III. PROPERTIES OF 
LINK-TRANSFORMATIONS 

A. Local operations 

Suppose that on one qubit, say qubit b for deflniteness, 
we perform a general local operation, not necessarily uni- 
tary Ub i.e. 



Mb = tla[{Ma(E)Ib){Ia<E)Ub) Pab (la 



(14) 
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In a cycle that includes the links a ^ b and b ^ c, this 
transformation would change both S{b,a) and S{c,b). 
For example, S{b, a)ji would be transformed into 



to find 



1 



(15) 



We can write this local operation on b as the left action 
on S{b, a), 



S{b,a) -> U{b)S{b,a) 



(16) 



where the components of the new matrix U{b) are given 
by 

[/(%,, = itr . (17) 

One can make a similar local operation, Ua, simultane- 
ously on qubit a and find that 



S{b,a) U{b)S{b,a)U{ay 



where 



(18) 



(19) 



Under local operations we see that the link- 
transformations change in the same way as the parallel 
transporters in lattice gauge theories if the gauge group 
is an orthogonal group. The total transformation around 
the closed loop described by eq. under arbitrary 
local operations, therefore becomes 

S{a,a]C)= (20) 
U{a)S{a, z)U{zf ■ ■ ■ U{c)S{c, b)U{bfU{b)S{b, a)U{af. 

Provided the local operations cancel each other, the trace 
of S{a, a;C) is invariant under these operations. We now 
prove this fact specifically for U e SU{2), that is U'^U = 
I. Each of the components of U'^U is given by 

UF,i,Ui,i, = itr {Uai.U^ai,) tr {U^ai,Uai,) . (21) 

Writing each of the 2x2 matrices in index notation we 
have 

-Ui^i^ {(Til (^'2 )MiM]ij2 i'^h )j2j-Mj-iji {^h )]4ji 

(22) 

where in the last equation the I indices take the integer 
values ... 3 and the i and j indices take the integer 
values and 1. Summation is implied by a repeated 
index. We can use the relation 



1 ^ 

9 X!('^')*1^2('^')»3»4 = ^lllji; 



(23) 



/=0 



(Si=o cTi = SWAP) and the unitary property of U 

U^^^MU^l4l^M^.^.=^n^. (24) 



tr i<7i^cTi^) = (5/^/3 



(25) 



The remaining local unitary transformations, those made 
at the beginning (or end) of the loop cancel from the 
cyclic property of the trace. So the quantity tr S{a, a; C) 
is indeed invariant under all local unitary transforma- 
tions. 

A simpler way to see that the local unitary transforma- 
tions do indeed cancel is to recognize that an arbitrary 
unitary acting on a qubit when written in terms of the 
Pauli matrices is simply a three dimensional spatial ro- 
tation acting on the three spatial components cri , (72 and 
CT3, that is, they are just rotations of the Bloch sphere. 
In other words the local operations U{a), U{b) acting on 
qubits a and b respectively in the 5(6, a) basis can be 
written explicitly as 



U{b)S(b, a)U{a) 



where and Ra are 3x3 rotation matrices, elements 
of S'0(3). That is, in the correlation matrix basis, 
U e SU{2) become U £ 5*0(3) due to the well known 
homomorphism SU{2) ~ SO{3) The components 

Sji are expectation values of the local spin measurements 
af and made on pab- One can verify the form of Ra 
and Rb using eqns. pT|) and 

In a similar way one can see that invariants where one 
performs a spin flip on each and every qubit are invari- 
ant under 14 e 5i(2,C) representing general local qubit 
operations up to a positive constant. The total trans- 
formation obtained by spin flipping every qubit can be 
explicitly written 




5(5,, a;C) = 5(a, z)r/ ■ ■ ■ r/S{c, b)riS{b, 0)77. 



(27) 



Under local operations U E 5i(2,C) we have seen from 
equations (|16m9|) that the correlation matrices 5 trans- 
form as 5(5, a) — >■ U{b)S(b,a)U{a)'^ thus we can form 
products such as U{a)rjU{a)^ from a transformation 
around a loop. Provided 



(28) 



our spin flipped quantities I{ab---z) are invariant. In 
fact eq. (|28p is the defining property of the group of 
Lorentz transformations, 50+(l,3), and due to the well 
known homomorphism 5L(2,C) ~ 50+(l,3) it indeed 
turns out that in the correlation matrix representation 
U e 5L(2,C) becomes U € 50+(l,3) [l|]. One can 
verify eq. (^5)) holds explicitly using eqs. PT)) and p^ . 
Therefore the spin flipped quantities tr5(a, a; C) are in- 
variant under local SL(2, C) operations. 
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B. Directional property 

One other useful property of the correlation matrices or 
link-transformations is simply demonstrated: The link- 
transformation taking m° to m^, the real matrix S{b, a), 
is the transpose of the link-transformation taking m'' to 
m°. That is 



S{b,a) = Sia,by 



(29) 



This property is easily seen from eq. dH). We note 
that this is another property shared by the parallel 
transporters in lattice gauge theory, the parallel trans- 
porter that takes you from one lattice point to another 
is the transpose of the parallel transporter that takes 
you back provided the gauge group is 0{N). How- 
ever, the parallel transporters have the additional fea- 
ture that a loop not enclosing area is the identity, for 
example [/(a, b)U{b, a) = 1. A similar expression for link- 
transformations does not hold. In fact we will make use 
of this property in the following section. 



IV. IDENTIFICATION OF INVARIANTS 

For any collection of qubits, one can consider the man- 
ifold representing the set of orbits of pure states under 
all local unitary transformations. That is, each point in 
the manifold corresponds to such an orbit. For a system 
of three qubits — we call them a, b, and c — it is known 
that the manifold of orbits is five dimensional (A 
quick but incomplete counting argument makes this re- 
sult plausible. The eight-dimensional space of pure states 
can be parameterized by fourteen real numbers, if we fix 
the normalization and the overall phase. A generic orbit 
has nine degrees of freedom, because each of the three 
local unitaries has three real parameters. This leaves five 
parameters to specify the orbit itself.) Stating this in an 
alternative way for the case of a pure three qubit state, 
\ip) is locally equivalent to {ip) provided all local invari- 
ants specifying the orbit are equal [l^. In the case of 
equality 1-0) and |0) have the same entanglement prop- 
erties and one can obtain |0) from \^p) simply by making 
local unitary transformations on each qubit. 

Several authors have studied local invariants of pure 
three qubit states [1-3, 7, 19-23]. In particular, Sudbery 
3 has identified a convenient set of algebraically inde- 
pendent invariants, each of which is a polynomial in the 
eight complex components of the vector !■(/;). Not count- 
ing the normalization (which is Sudbery's /i), there are 
five of these invariants, the same as the number of di- 
mensions: 

/2-tr[(p,)2], /3=tr[(pfc)2], /4=tr[(p„)2] 

— {^Pab)ij' ,i' j {Pbcjjk'.j'k {^P ca) ki' ,k' i 

h - (l/4)rf,, (30) 

where summation over repeated indices is implied in the 
definition of /s. Here each index takes the values and 



1, and we have used the letters i, j, and k to refer 
to qubits a, b, and c respectively. The Kempe invari- 
ant /s 4] can be written in several different ways, the 
above form being most convenient for our purpose. The 
quantity Tate, is the 3-tangle which measures a kind of 
three-way entanglement characteristic of the GHZ state 
(1/a/2)(|000) + |111)) [11. If we write in terms of 
the standard basis states as \'ip) = J^'^ijkl'ijk), then the 
invariant /g can be expressed as 

(31) 

where eij is the antisymmetric tensor in two dimensions. 

The invariants listed in eq. (j30p are not complete in 
the sense of determining a unique orbit. In particular, 
these invariants do not distinguish between a state and 
its complex conjugate, which may well lie on different 
orbits. Because Ii - Iq are real 



(32) 



As reported by Acin et al. [2J|, Grassl has shown that 
this ambiguity can be removed by including a single bi- 
nary invariant based on a complex twelfth-degree poly- 
nomial in the amplitudes aijk- 

We now ask whether the invariants in eq. ([50)1 can be 
generated via the formalism of section [Hi The first three 
can indeed be expressed quite simply in this way. For 
example. 



liab) = ti-{S{a,b)S{b,a)} 




(33) 



The last line follows from the fact that the operators 
{af (8) o'°)/2 constitute a complete orthonormal basis for 
the space of 4 x 4 matrices. 

The Kempe invariant I5 fits particularly well into our 
scheme. As we now show, this invariant is simply 



I5 = I{abc) = tr {S{a, c)S{c, b)S{b, a)} . 



(34) 



To see this, we start with the following expression for 
I{abc): 



1 



tr [{^i ® ^l)Pab] tr [{al ® a^pbc] tr [« ® a^p. 



(35) 

in which summation over Z, to, and n is implied. This 
summation considerably simplifies the expression, be- 
cause of eq. This relation tells us how to connect up 
the indices of the three density matrices, and we obtain 
the interlocking pattern that we saw in eq. (pO)) : 



I{abc) = {pab)ij',i'j iPbc)jk',j'k [Pcajki' ,k' • 



(36) 



Of the set of invariants that Sudbery identifies, the 
only one remaining is /g, a degree-eight polynomial in 
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(a)/(afe) 



o o 



o 

{h)I(bc) 



gradients at any point, together with the gradient of the 
normaUzation invariant Ii = aijkO-*jki constitute a Un- 
early independent set of vectors [3|. One finds that this 
is indeed the case. So we now have the following list of 
path-based invariants, not quite identical to Sudbery's 
but no less complete: 



(c)/(ca) 



I{ab) = tv{S{a,b)S{b,a)} , 

I{hc) ^ tr{S{h,c)S{c,h)} , 

I{ca) — ir {S{c,a)S{a,cy\ , 

I{ahc) = tr{S{a,c)S{c,h)S{h,a)} , 

I{abab) = ir{S{a,b)S{b,a)S{a,b)S{b,a)} , 



(40) 



o • 

(A)I{abc) (e)I{abab) 



FIG. 2; Graphical representations of the closed paths giving 
the set of polynomial local SU (2) invariants for pure states of 
three qubits. 

the components Uijk and their conjugates which is pro- 
portional to the square of the 3-tangle. Iq is also invari- 
ant under SL(2, C) unlike Ii - which do not have this 
higher invariance. Our formalism does not produce /g 
directly, though we can easily generate a different poly- 
nomial of degree eight that is likewise algebraically inde- 
pendent of the first four invariants. It is defined by any 
path that cycles twice between two of the qubits, that is, 
any of the invariants 

I{abab) = I{bcbc) = I{caca). (37) 

One can write down this invariant directly in terms of the 
reduced density matrices pab and pba (which are related 
to each other by the swap operation that interchanges the 
two qubits), following precisely the pattern of index con- 
nections that we see in eq. (j36p . Now, however, because 
the same two-step path is repeated, we use subscripts 1 
and 2 on the indices to distinguish the two round-trips. 

I{abab) = tr {S{a, b)S{b, a)S{a, b)S{b, a)} 

= iPab)iij[,i[ji {Pba)jii'^,j[i2 iPab)i2j^,i'2j2 (Pba) j2i{,j^i!^) 

We can alternatively write out this invariant in terms of 
the components aijk'. 

I{abab) = atij2kia*23ikiO'U3ik2^*i3j2k2 

In this latter form it is clear that the invariant is sym- 
metric under permutations of the qubits: by permuting 
the factors of a and a*, one can interchange the roles of 
the i, j, and k indices. 

To show that the invariants we have identified are al- 
gebraically independent, it is sufficient to show that their 



the last one being symmetric under permutations of the 
three qubits even though the path it is based on is not. 
The three kinds of path we have used in constructing our 
invariants are illustrated in figure [H 

Notice that our construction provides an interpretation 
of the Kempe invariant I{abc) = I^- In our measurement 
scenario, in which each successive measurement collapses 
the state of the next qubit, the Kempe invariant is the 
trace of the transformation that results from following 
the triangular path through all three qubits. Recall that 
at each stage in this measurement scenario, the trace of 
the new M matrix is equal to the probability of getting 
the desired outcome. Thus, the resulting invariant tends 
to be larger if the collapsed state at each step is strongly 
represented in the original reduced density matrix of the 
qubit in question. The most extreme example of this 
kind of consistency is the case of a completely factorable 
state, in which the collapsed state must be proportional 
to the original pure state of the given qubit. And indeed, 
the Kempe invariant is largest when the state is fully 
factorable {I{abc) = 1). One can also show the Kempe 
invariant takes its minimal value for the W state |00I) -I- 
|010) + |100) at I{abc) = 2/9 lOj. 



V. LOCAL INVARIANTS USING THE SPIN 
FLIP 

In section IIVI we provided a complete set of alge- 
braically independent local SU{2) invariants by consid- 
ering different closed paths around the lattice. In this 
section we again consider closed paths around the lattice 
but this time including the spin fiip operation on every 
qubit lifting the invariance of the quantities produced to 
5L(2,C). 

From the SU{2) invariant list, it turns out that we 
can replace I{ab), I {be), and I{ca) with I{ab), I {be), and 
I{ca); the invariants are still independent. Moreover, 
these 'flipped' invariants can also be interpreted in terms 
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of entanglement. One finds that 

I{ab) = tr [pabPab] = Tab + ^Tabc 
I (be) = tr [pbcPbc] = Tbc + -jTabc 
I(ca) = tr [pcaPca] = Tea + ^Tabc (41) 

Here p = [a2® (^2)9^ {'^2 ^ (^2) is the spin flipped state 
of the two qubits, and Tab is the tangle between qubits a 
and b, a measure of their pairwise entanglement (it is the 
square of the concurrence) 25] . The completely flipped 
version of I{abab), that is, I{abab), likewise produces a 
non-trivial invariant but it is not algebraically indepen- 
dent of I{ab), I (be), and I{ca). A natural eighth order 
SL(2,<C) invariant is given by the determinant of any of 
the link-transformations, for example: 

det [S{a, b)] = -^Tab{Tab + Tabc)- (42) 

16 

One can see the determinant of the link-transformations 
is indeed SL{2,C) invariant from the way the link- 
transformations change under arbitrary local SL{2, C) 
transformations (eq. (jlSp ) and the property of the de- 
terminant det[AB) — det (A) det (-B). We can show this 
invariant is algebraically independent of /(a6), I {be), and 
I{ca) using the same methods of section ITVl 

So far we have a set of four SL{2, C) invariants. They 
tell us about the entanglements in the state since one 
can reconstruct the amounts of entanglement. Tab, t^c, 
Tac and Tabc from just these four invariants. For example, 
the 3-tangle can be expressed as 

Tabc = 2y'l6 det S'(a, fe) + /(a6)2. (43) 

Even though the three qubit labels do not enter this ex- 
pression symmetrically (there is no explicit reference to 
qubit c), the 3-tangle is symmetric under permutations 
of the qubits. Similar expressions can be written for the 
2-tangles. 

These four amounts of entanglement do not form a 
complete set of algebraically independent invariants. To 
complete this set we could use the SU (2) Kempe invari- 
ant which one can verify is algebraically independent of 
the four tangles, however we would also like an invari- 
ant with SL{2,C) local invariance. A natural choice of 
loop is the one that gave the Kempe invariant. Strangely 
enough, the completely flipped version of the Kempe in- 
variant, that is, I{abc), turns out to be exactly zero for 
all pure three qubit states as we prove in the following 
subsection. 




{d)I{abc) (e)det(5(a, fe)) 



FIG. 3; Graphical representations of the closed paths giving 
the set of polynomial local SL{2, C) invariants for pure states 
of three qubits. Notice that I{abc), the flipped Kempe invari- 
ant, is the only closed path to enclose area. It is also always 
equal to zero. 

A. Proof of I{abc) — for any pure three qubit 
state. 

From the definition of I{abc), one finds directly that 

I{&k) =1 - tiipl) - tr(p2) _ tr{pl) 

+ tr[{pa (g) Pb)Pab] + tr[(pt (g) Pc)pbc] (44) 

+ tr[(pa (g Pc)Pac] - h- 

Sudbery Q showed that Kempe 's invariant can also be 
written as 

h = 3tr[(pa ® Pb)pab] - tr(p3) - tr{pl) 

= 3tr[(p6 ® pc)pbc] - tr(pg) ~ tr(p3) (45) 

= 3tr[(pa ® pc)pac] - tr(p3) _ tr(p3). 

Using these relations we can rewrite I{abc) as 

I{ab~c) =1 - MpD + tr(pg) + tr{pl)] 

2 , . . (46) 

This last expression is a function only of the trace of 
powers of the single-qubit density matrices. The Cayley- 
Hamilton theorem for any 2x2 matrix X is 

A2-tr(A)A + dct(A)I = 0. (47) 

Multiplying this expression by X, taking the trace, and 
then using the fact that for a single qubit, det p — 
(1/2)(1 — trp^), we obtain the relation 

l/3-tr(p2) + (2/3)tr(p^) =0, (48) 

which together with eq. (^5)) shows that I{abc) — 0. 
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VI. OPERATIONAL INTERPRETATION OF 
INVARIANTS 

Going back to our original thought experiment in sec- 
tion HI] we can obtain a rigorous operational interpreta- 
tion of our invariants as follows if we drop the restriction 
that Ma is positive but still Hermitian. In other words 
Ma represents an observable rather than a measurement 
outcome. If we relax the positivity our invariant can be 
thought of as the average fidelity between our initial ob- 
servable outcome Ma and our final observable outcome 
M'^ resulting from the transform around the loop asso- 
ciated to the particular invariant, /(C). The average is 
taken over all possible initial outcomes Ma with a fixed 
size. Making this idea more precise, we have 

/(C)cx(tr(MtA<)), (49) 

where the brackets (.) denote the average. 
The constraint on Af^'s size is given by 

tv{MlMa) = 2fc2. (50) 

In terms of the Pauli operator basis, this condition can 
be written as 

3 

tr{MlMa)=2Y,{<f ^'^k\ (51) 

1=0 

The must be real for Ala to be Hermitian and rep- 
resent the outcome of an observable. Similarly we can 
write the outcome on a following the loop C in the Pauli 
operator basis as (dropping the sub- and superscript a) 

3 

M' ^^m[a, (52) 

i=0 

where 

m[ = S{C)ijmj. (53) 

S{C) is the total transformation around the loop. 

We can substitute these expressions into the equation 
for the fidelity between transformed and initial observ- 
able outcomes to give 

3 

tr(M^M') 

1=0 

3 

= 2 miS{C)ijmj. (54) 
i,i=o 

Our invariant /(C) is given by the elements S{C)u and 
therefore we want to find an expression solely in terms of 
these elements. 

We now average this fidelity. Since eq. (|50p is the equa- 
tion of a 3-sphere with radius k we can perform the av- 
erage over the surface of the 3-sphere. Writing the in 



hyper-spherical coordinates we have 

mo = k cos (j)i 

mi = k sin (pi cos 02 

m2 = k sin (j)i sin 02 cos 03 

m3 = fc sin 01 sin 02 sin 03 . (55) 

We can now compute the average fidelity in terms of these 
coordinates. It is given by 

3 

(tr(A/^A/')) ^"^Y] S{C)ij [ rmmjdS. (56) 
^ ^,J=o ^5 

S is the entire surface of the 3-sphere, dS = 
A;^ sin^ 01 sin 02^01 (i02d03 is the area element and A — 
27T^k^ is its total surface area. One finds that 

4 / m.nijdS = y(5y (57) 

and we obtain our desired result 

(tr(A/tA/')) = y/(C). (58) 

One can also average over all observable sizes (or 
strengths) k to obtain the same result up to a constant. 
The surface integral over the sphere now becomes a vol- 
ume integral over the 3-ball. We also note we can choose 
Ma to be an element of SU{2). That is, a unitary that 
does not have to be Hermitian. In this case k = 1 and 
the average is over the three Euler angles describing a 
element of this group. For SU(2) mg is real and mi, m2, 
ni3 are purely imaginary. The proof goes through in the 
same way. 

VII. CONCLUSIONS 

In this paper we have presented a geometric approach 
to constructing quantities that are invariant under lo- 
cal SU{2) and SL{2,C) transformations. Our basic con- 
struction corresponds to a scenario in which a measure- 
ment outcome on each particle along a closed path de- 
fines the state of the next particle. We have seen that 
one can produce in this way an algebraically indepen- 
dent set of five invariants for a pure state of three qubits, 
almost identical to the set of invariants identified by Sud- 
bery. One of these quantities, the Kempe invariant, has 
been difficult to interpret as an amount of entanglement. 
In our construction, though, it is the one that emerges 
the most naturally. Unlike the other four invariants, the 
Kempe invariant = I{abc) corresponds to a path that 
'encloses area' in the sense that one does not retrace one's 
steps. This property sets the Kempe invariant apart from 
the others. Notice that for an area enclosing path, one 
needs at least three qubits. In a future paper we will 
exploit this area enclosing property and the existence of 
a special form of a polar decomposition for correlation 
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matrices to find quantities much more analogous to lat- 
tice gauge field theories. The gauge group will turn out 
to be the group of Lorentz transformations and has an 
operational interpretation in terms of general local qubit 
operations. The invariants, the Wilson loops, in this con- 
struction, will be related to the curvature of the correla- 
tion space [25|. 

We have also provided an operational interpretation of 
the invariants, including the Kempe invariant, in terms 
of the average fidelity between initial and transformed 
observable outcomes. 

We have concentrated on pure three qubit states as a 
test ground for our ideas, however there is nothing spe- 
cific here about the numbers of qubits of our quantum 
state. Neither is there any requirement for the state to 
be pure or for the subsystems to be two level. The frame- 
work presented here can be applied to any qudit state 
to generate local unitary invariants. To make this gen- 
eralization, one replaces the Pauli matrices specific for 
qubits, by the generalized Gell-Mann matrices Ai, an or- 
thonormal basis for the (real) {(P — 1) dimensional vector 
space of traceless hermitian dx d matrices with the inner 
product {X,Y) = tT{XY). For U e SU{d), the com- 
ponents of the local operation in the correlation matrix 
basis, Uij = triUXiU^ Xj) is a special orthogonal matrix 
in SO{d^ — 1), representing U in the (d^ — 1) dimensional 
(adjoint) representation of SU{d). More precisely, the or- 
thogonal matrices with components Uij form a subgroup 
of SO{d^ ~ l)i being the homomorphic image of SU(d) 
called the adjoint group, which is the quotient of SU(d) 
by its centre, the subgroup of matrices ojld where uj is 
a dth root of unity. Thus, the property U^U = Id still 
holds and tr5'(C) is invariant under SU{d). However, for 
SL{d, C) the connection with a Lorentz group only works 
for d = 2. 



[1] H. A. Carteret, A. Higuchi, and A. Sudbery, J. Math. 

Phys. 41, 7932 (2000). 
[2] N. Linden and S. Popescu, Fortschr. Phys. 46, 567 

(1998). 

[3] A. Sudbery, J. Phys. A: Math. Gen. 34, 643 (2001). 
[4] J. Kempe, Phys. Rev. A 60, 910 (1999). 
[5] M. Grassl, M. Rotteler, and T. Beth, Phys. Rev. A 58, 
1833 (1998). 

[6] E. M. Rains, IEEE Trans. Inf. Theory 46, 54 (2000). 
[7] H. Barnum and N. Linden, J. Phys. A: Math. Gen. 34, 
6787 (2001). 

[8] M. Teodorescu-Frumosu and G. Jaeger, Phys. Rev. A 67, 
052305 (2003). 

[9] A. Osterloh and J. Siewert, Phys. Rev. A 72, 012337 
(2005). 

[10] A. Osterloh, Appl. Phys. B 98, 609 (2010). 
[11] G. Miinster and M. Walzl, hep-lat/0012005 (2000). 
[12] R. Mosseri and R. Dandoloff, J. Phys. A: Math. Gen. 34, 
10243 (2001). 

[13] B. A. Bernevig and H.-D. Chen, J. Phys. A: Math. Gen. 
36, 8325 (2003). 



We have seen for three pure qubits that only one area 
enclosing path exists, but as the number of qubits in- 
creases, there should be many more Kempe like quan- 
tities, since there will be many more paths that enclose 
area. Our original construction, on the other hand, pro- 
duces many more invariants even for three qubits, be- 
cause it makes critical use of the 'shrinking' component 
of the transformation defined by the spin correlation ma- 
trix. 

One caveat with this approach as it stands is that we 
have only used the information contained in the two qubit 
density matrices. For states with large numbers of qubits 
one cannot obtain a full set of invariants since too much 
information about the overall state is lost when tracing 
out all but the two qubits in each link. This approach 
could be extended by considering contractions of the full 
correlation tensor. For example, an N qubit system is 
described by a real tensor Si^j^-.-z^ where «i, j2, • • • ,zn 
take the values 0,1,2,3. It would be interesting to find 
an operational interpretation of contractions of these 
A^— tensors. 

Acknowledgments 

We thank Johan Aberg, Stephen Brierley, Caslav 
Brukner, Berge Englert, Richard Jozsa, Noah Linden, 
Ognyan Oreshkov, Jiannis Pachos, Wonmin Son and An- 
dreas Winter for helpful comments and discussions. For 
financial support MSW acknowledges EPSRC, QIP IRC 
www.qipirc.org (GR/S82176/01), the NRF and the MoE 
(Singapore) and an Erwin Schrodinger JRF. ES and VV 
acknowledge the National Research Foundation and the 
Ministry of Education (Singapore). ME acknowledges 
support from the Swedish Research Council (VR). 



[14] P. Levay, J. Phys. A: Math. Gen. 37, 1821 (2004). 

[15] P. Levay, Phys. Rev. A 71, 012334 (2005). 

[16] A. Jamiolkowski, Rep. Math. Phys. 3, 275 (1972). 

[17] F. Verstraete, J. Dehaene, and B. D. Moor, Phys. Rev. 

A 64, OIOIOI(R) (2001). 
[18] P. Arrighi and C. Patricot, J. Phys. A: Math. Gen. 36, 

L287 (2003). 

[19] V. Coffman, J. Kundu, and W. K. Wootters, Phys. Rev. 

A 61, 052306 (2000). 
[20] A. Acin, A. Andrianov, L. Costa, E. Jane, J. I. Latorre, 

and R. Tarrach, Phys. Rev. Lett. 85, 1560 (2000). 
[21] R. M. Gingrich, Phys. Rev. A 65, 052302 (2002). 
[22] M. S. Leifer, N. Linden, and A. Winter, Phys. Rev. A 

69, 052304 (2004). 
[23] T. A. Brun and O. Cohen, Phys. Lett. A 281, 88 (2001). 
[24] A. Acin, A. Andrianov, E. Jane, and R. Tarrach, J. Phys. 

A: Math. Gen. 34, 6725 (2001). 
[25] W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 
[26] M. S. WiUiamson, M. Ericsson, M. Johansson, 

E. Sjoqvist, A. Sudbery, and V. Vedral, arXiv:1102.5609 

(2011). 



10 



[27] For the case of pure three qubit states wc need five contin- mial invariant to identify which states are locally equiv- 

uous local polynomial invariants and one binary polyno- alent. This will be discussed later in the section. 



